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Abstract 



A Classical Trajectory Monte Carlo (CTMC) simulation has been made 
of processes of charge exchange and ionization between an hydrogen atom 
and fully stripped ions embedded in very strong static electric fields (O(10^° 
V/m)), which are thought to exist in cosmic and laser-produced plasmas. 
Calculations show that the presence of the field affects absolute values of the 
cross sections, enhancing ionization and reducing charge exchange. Moreover, 
the overall effect depends upon the relative orientation between the field and 
the nuclear motion. Other features of a null-field situation, such as scaling 
laws, are revisited. 

To be published in Journal of Physics B 
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I. Introduction 

The study of the effects of externally imposed electric or magnetic fields 
on atomic systems has been since a long time an active area of research in 
physics, both when dealing with isolated atoms (Stark and Zeeman effects, 
photoionization, chaotic spectra, modifications of energy levels and wave, see 
e.g. Seipp and Taylor 1994, Wang and Greene 1994, Delande and Gay 1986, 
Priedrich 1990, Zang and Rustgi 1994, Delande et al 1994, Bylicki et al 1994, 
Buchleitner et al 1994), as when the attention is focused onto many-particle 
processes (electron-atom or atom-atom collisions). 

In order to influence the electronic dynamics it is necessary that the force 
exerted by the field is comparable with the Coulomb intra-atomic forces, 
which are quite strong if the atom is not in a highly excited state: for mag- 
netic and electric fields the natural units of measure are Bq = mle^c/h^ = 
2.35 • 10^ T, and Fq = mle> /h^ = 5.14 ■ 10" V/m. Fields of such strengths 
may be generated only under rather exotic situations: in white dwarfs mag- 
netic fields of 10^ T, and in neutron stars up to 10^ T, have been discovered 
(Ruder et al 1994). The huge unipolar induction currents which generate 
the magnetic field are sustained by differences of potential of order 0(10^^ 
V/m) (Ruder et al 1994). By a different mechanism, electric fields may be 
generated even independently by magnetic fields: it is the double layer (DL) 
phenomenon which is expected to be found, for example, on the surface of 
neutron stars (Williams et al 1986, Raadu 1989). A DL is a local discon- 
tinuity surface consisting microscopically of two oppositely charged layers 
of plasma. The resultant electric field within this region is very strong and 
acts to re-establish a neutrality condition which is violated, for example, as a 
consequence of temperature gradients on the surface of the star or because of 
charge density gradients produced by the different rates at which electrons 
and ions, falling on the star surface from the outer space, are decelerated 
(Williams et al 1986). It is thought that fields > 10^° V/m may be generated 
through this mechanism. 

DLs of considerable strength are created also in Earth laboratories during 
laser experiments devoted, for example, to inertial confinement fusion studies: 
in these experiments a small target of hydrogen isotopes is irradiated with 
high power 10^^ ^ 10^° W/m^) laser pulses. The sudden heating, which 
converts into a hot plasma most of the target, drives a fast dilatation of the 
gaseous outer corona where, because of their higher velocity, the electrons 
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will tend to lead the expansion. The result is a DL which slows down the 
electrons and accelerates the ions (for a review of the subject see Eliezer and 
Hora 1989). Numerical simulations show that, inside these DLs which have 
a width of about 10~^ 10~^ m, electric fields of 10^ lO-*^^ V/m may be 
reached (Eliezer and Hora 1989, Eliezer et al 1988). 

Charged particles within a DL usually have great kinetic energies because 
of the heating, the electrostatic force, and — in stars — the gravitational at- 
traction: they may be O(10KeV/amu ) (Williams et al 1986, Eliezer and 
Hora 1989). DLs thus provide us with rather extreme conditions where to 
examine interatomic processes, particularly collisions. In this paper we aim 
to perform a numerical investigation of the processes of electron capture and 
ionization between a neutral atom and a bare ion. This kind of study is 
not new in literature: a similar work has been performed by Grosdanov and 
McDowell (1985) and McDowell and Zardona (1985) for very strong, astro- 
physical magnetic fields. Electric fields have been instead studied by Olson 
and MacKellar (1981). There, however, the attention was focused on fields 
attainable in ordinary laboratory conditions, i.e. much weaker than those 
considered here; to compensate for the smallness of the field only Rydberg 
atoms had to be considered. Here we carry on the analysis on low lying 
electronic states instead that highly excited ones; some topics not covered 
by Olson and MacKellar (1981) are touched: the results of the scattering 
process, expressed in terms of cross sections, are studied with respect to the 
charge of the incident particle and the initial state of the system. Also, par- 
ticular attention is paid to the geometrical features (direction and spatial 
extension) of the external field. 

II. Theory 

The collision process has been studied using the CTMC method, in which 
both nuclei and the electron are considered classical particles obeying Newton 
laws but the initial conditions of the electron are randomly chosen within a 
distribution which simulates the quantum mechanical behavior. The study of 
the process is then reduced to numerically integrate a set of coupled ordinary 
differential equations. 

The CTMC method, developed by Abrines and Percival (1966a), has been 
frequently used in studies of atomic scattering, being quite easy to implement 
and rather accurate in results. The method works at its best in the region 
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of intermediate impact energies, i.e. when the relative impact internuclear 
speed Vp is greater than the classical orbital electron velocity Ve- Vp > Ve] 
in this work we will consider the range 1 < Vp/v^ < 3.5 (corresponding to 
energies E < 300 KeV/amu), which is quite in accordance with the velocities 
expected to be found in DLs. 

The equations of motion for each of the particles are derived from the 
Hamiltonian 

3 T)2 3 
^=E|+ E (1) 

1=1 ^ i<j=i iy« yji i=i 

where Qi, Pi {i — 1, 2, 3) are the coordinates and the momenta conjugate for 

— * 

the two nuclei and the electron; Zi the charge of the i-th particle, and F the 
electric field (here and in the following, for the sake of easiness, the electric 
field will be expressed in atomic units; in these units Fq of section one is 
equal to 1). 

An atom embedded in a strong electric field is an instable system: classi- 
cally, a bound electron is field-ionized if the strength F of the field is greater 
than E^/i, with E^ binding energy. Thus, any classical calculation loses 
meanings at strengths beyond this value (for an hydrogen atom in the ground 
state, E^/i = Fq/16 = 0.0625 au). However, quantum mechanically, there is 
not a ionization threshold, instead all bound states turn into resonances of 
finite lifetime, whatever small F be. Several calculations of lifetimes exist in 
literature; here we shall refer to Damburg and Kolosov (1976): in that work 
it is estimated that, in a field of F = 0.04 au the mean life of an electron in 
the ground state is ~ 2- 10~^^ s. A particle moving at a speed Vp > 10^ m/s 
(corresponding to a kinetic energy > 25 KeV/amu) travels during this time 
a distance Ip — VpTg > 10~^ m. If we refer, for example, to a laser-produced 
DL, its characteristic spatial length is given by the Dcbyc length Xd; usually, 
this value is < 10~^ m (sec Eliezer and Hora 1989, sect. 4). 

Further, in a laser-produced DL the electric field is not at all static, the 
timescale needed to estabhsh it and over which its variations are not negligible 
being in the range of 10~^^-=-10~^^ s (Ehezer and Hora 1989, sects. 4,5). Even 
within such short times, a particle at the speed Vp travels for several hundreds 
of A, which is well beyond the typical length over which charge transfer 
processes take place. So we see that, even under rather extreme conditions, 
a classical description of the coUision process does not lose its meaning (we 
remark however that for cosmic DLs the picture is not so satisfactory, because 
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the width of the DL is supposed to be greater, thus in this case the vahdity 
of our hypotheses may be questionable). 

Two further remarks need to be made when deahng with excited atoms 
(we will consider only n — 2): first, the validity of the classical approximation 
still holds because we shall consider weaker fields, ad hoc rescaled to com- 
pensate for the smaller binding energy. Further, because of the Stark mixing 
between levels due to F, and in particular between the levels 2s — 2po, all 
the electrons are allowed to decay by spontaneous emission of radiation into 
the ground state within short times (~ 10~^ s) which, however, are always 
longer that the timescales of the collision process. 

The electric field exerts also an influence on electronic energies (Stark 
shift) and wave functions. As far as the first are concerned, we see from 
Damburg and Kolosov (1976) that the correction to the null- field values is 
about 1 percent and thus may be neglected in the following. 

Using a bit more care instead is needed when dealing with position and 
momentum distributions; in CTMC method quantum wave functions are 
replaced with a random sampling of the initial coordinates from statistical 
distributions. The older and more used choice is the microcanonical distri- 
bution (Abrines and Percival 1966b), where momentum pe and position 
are picked up from an uniform distribution subject to the bound that the 
total energy is equal to its quantum mechanical value. 

This method provides a statistical momentum distribution which is equal 
to the quantum mechanical one, while the radial distribution is rather bad, 
showing a cut-off beyond r = \qe\ = 2n^ au . This causes a wrong estimate of 
total cross section at low impact energies, underestimating the contribution 
from coUisions at large impact parameters. Several methods have been de- 
vised to overwhelm this deficiency, with good results (see Hardie and Olson 
1983, Cohen 1985, and also Eichcnauer et al 1981, Kunc 1988, Reinhold and 
Falcon 1988a,b, Schmidt et al 1990 for a general discussion over the CTMC 
and other classical approximations). However, any classical method which 
attempts to simulate a correct radial quantum distribution has to give up 
the assumption that electron energy is a well fixed quantity. The electron 
will thus have some probability to be given an energy beyond the ionization 
threshold when the field is turned on, and to be field ionized. This spurious 
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contribution will add to the correct probability of ionization for ionic im- 
pact. As an example, in Cohen's model (Cohen 1985), with a field strength 
of Fo/2, about 30 per cent of the electrons would have an energy beyond the 
ionization threshold because of the choice of the initial conditions. 

Initial position and momentum have been chosen as if the electron were 
not under the influence of the field. In Fig. 1 wc report the evolution of 
spatial distribution; it is clear that already after a short period, < 5 au, the 
spatial distribution has reached the equilibrium value compatible with the 
field. 

Target and projectile nuclei have been set at an initial distance of 20Z 
au, where Z is the projectile charge. The equations of motion have been 
integrated until the nuclei were well far apart. For each test a number of 
runs rangmg between 3-10^ and 26-10^ was performed. We did not attempt 
to minimize the statistical error for every point, however for all of the data 
the error goes from 1 to a maximum of 20 per cent. 

When considering the final state of the electron it is essential to remem- 
ber that an electron captured in an high enough quantum number may be 
subsequently field ionized. We have followed the convention of Olson and 
MacKellar (1981) of considering as ionized such electrons. 

Three different field geometries have been considered: adopting the sys- 
tem of reference where the neutral atom is initially at rest and the ion is 
moving along the z axis towards the positive direction, the electric field has 
been chosen either parallel to z, antiparallel or perpendicular. In the follow- 
ing we will adopt the convention of designate a parallel or antiparallel field 

— * 

according to its projection on the z axis, so F > means F \\ Vp. 
III. Hydrogen— proton collisions 

The first series of calculations was about H — collisions: in Fig. 2 we 
plotted electron capture cross section, acxi and ionization, cxjo^, versus impact 
velocity for two different values of the field, chosen parallel. It may be stated, 
alike shown in Olson and MacKellar (1981), that F enhances ionization, but 
at the same time opposes to the process of recombination on the projectile, 
working against charge exchange. The overall effect is quite clearly increasing 
with the strength of F as far as (Ticm is concerned, while for a^x the effect is 
less marked, even though visible. We note further that the field is effective 
on electron capture only at low energies, its infiuence being negligible above 
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f = 2 au (i.e. E = 100 KcV/amu). Instead, Uion is affected in tlie same 
manner along all the energy range. 

In Fig. 3 the same calculations have been performed by varying the 
field direction. We may state that is dependent from the sign of F, the 
difference being small but regular. In Olson and MacKellar (1981) it was 
already suggested that when the field pulled the electron in the direction of 
the outgoing projectile, electron capture and loss were augmented. A less 
definite effect is visible on aion- There are not appreciable differences when 
the field is set orthogonal to the ion direction or parallel to it. 

The calculations were then repeated with a change in the initial condition, 
the electron being now set in the n = 2 state. From Fig. 4, where we have 
plotted the data corresponding to two opposite values of the field, one argues 
that the difference between acx{F < 0) and acx{F > 0) is clearly enhanced, 
while not a clear behavior along the entire energy range appears for (7jo„. 

It is well known that the, in field-free conditions, the collision process 
has certain symmetries, basing upon which some scaling laws can be inferred 
(see, for example, Janev 1991): all collisions between an hydrogen atom in 
the state n and an ion of charge Z lie on the same curve when one plots cross 
sections in terms of the reduced quantities 

a(E) = ^-^, E=^. (3) 

^ ^ n^Z' y/z ^ ^ 

Under non zero-field conditions the same relations should hold, provided the 
external field is scaled according 

In Fig. 5 we have plotted the same quantities of Figs. 3 and 4 in terms of 
these new reduced variables; both for acx and aion a single trend is discernible, 
but the spreading of the points is not negligible. 

As far as differential cross sections are concerned, we have plotted in Fig. 
6 the quantities da/db for two scatterings at the same impact velocity and 
opposite values of F. The effect of F, in both cases, is to reduce the effective 
range within which capture takes place, while ionization is now made possible 
at larger impact parameters. 

IV. Hydrogen— multiply charged ions scattering 
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In order to study features of non-resonant scattering we made simulations 
of collisions of hydrogen with He^^, Li^^, C®+, and O^^ ions, whose results 
are reported in Fig. 7. We may observe that a different behavior exist 
between the lightest ions (H, He), and the others. This is also clearly revealed 
in Fig. 8 where the same quantities have been plotted in terms of the reduced 
variables of eq. 3; however we observe that Z-scaling of eq. 3 still holds for 
the heaviest ions and on the whole this scaling law is obeyed quite well. 

A further subject is the distribution of acx over quantum states of the 
projectile. In this work we only considered distributions over the principal 
quantum number n. Classically, a bound electron may have any value of the 
energy, below the ionization potential. Capture into discrete energy levels is 
approximated by a procedure developed, for example, in Olson (1981) and 
Salop (1979). It is well known that there exist a proportionality relation 
between the ion charge Z and the quantum state n into which the electron 
is preferentially captured, n fa Z^^* (Olson 1981). In Fig. 9 we plotted these 
distributions for C®+ and O^^ ions: F has no effect when capture takes place 
in low lying states while, as one could foresee, capture in higher states, which 
are less bound to the nucleus, is strongly suppressed. The position of the 
maximum is unaltered with respect to the null-field case (Olson 1981, Salop 
1979). We remark that n = 7 and n = 9 are, respectively for carbon and 
oxygen, the maximum non field-ionized quantum numbers at this value of 
F. 

V. Shrinking the extension of the field 

Until now all the calculations have been performed under the hypothesis 
of an uniform external field, extending in all directions up to infinity. As a 

final issue we modified this configuration by limiting the field within a finite 
interval: it is abruptly set to zero outside a distance L from the target nucleus 
along the ^— axis in both directions. After this modification we have a more 
symmetrical situation, in which both Coulomb interparticle forces and F 
can exercise an influence over only a limited region. Ideally, the length L 
should be taken great enough to include most of the interparticle interaction 
within it: we chose L = 10 Z an . Repeating previous calculations for H—H^ 
collisions gives some modifications: is no longer uniformly suppressed and 
may become even greater than in the zero field case, and the opposite is true 
for a ion (Fig. 10). Another non neghgible difference stems when dealing with 
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high charged ions: in Fig. 11 wc have plotted n-distributions for H — C^^ 
impacts at a velocity Vp = 2. The absolute values of acx corresponding to 
F = 0, +0.03, -0.03 au are, respectively, 9.5 • 10"^^ cm^, 5.4 • lO'^^ cm^, 
17.4 • 10~^^ cm^, and the auxiliary contribution to the F < case comes 
from electrons which are bound in high lying states around n = 10, 11. We 
suggest that this result may be explained by a mechanism similar to capture- 
via-ionization (see McCartney and Crothers 1994): the incoming projectile 
strikes the electron, raising it in the continuum with an initial velocity vq, 
which we will assume to be parallel to the ion velocity Vp. The electron is 
then accelerated by the field. On z = L, it reaches the maximum speed 
which, with the position F = —\F\, may be written Vf ^ + 2L\F\. The 
relative distance between the ion and the electron is now 

D = v;^-L. (5) 

Capture takes place with preference when Vp and Vf closely match, Vp ~ Vf. 
The binding energy is then 



from which one gets n ^ [J ZD/2] ([ ] means the integer part). The condition 



Vp = Vf gives Vq = — 2L\F\. The substitutions Vp = 2, Z = 6, L = 
60, |F| = 0.03 give n 10 in good accordance with Fig. 11. A similar, but 
much weaker, feature appears in the F > case; the maximum at n = 7 may 
be partially explained by the same mechanism, provided that now the field 
has a braking effect on the electron, whose initial speed must thus be greater 
than Vp. 

VI. Summary 

In this work we have performed some numerical simulations with the 
aim to investigate how a strong static electric field may modify an energetic 
collision between an ion and an atom. We studied total and partial cross 

sections for charge-exchange and ionization versus some meaningful param- 
eters as impact velocity Vp, ion charge Z, initial target quantum state n and 
spatial extension of the field. It has been found that the field deeply affects 
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the results of the scattering; detailed effects depend upon the strength of the 
field and its orientation with respect the relative nuclear motion. Scaling 
laws, true in null-field situations, seem now more doubtful, and are surely 
true only within strict ranges of the parameters. 
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Figure Captions 



Figure 1: upper, radial distribution; lower, polar angle distribution. Full 
line represents micro canonical distribution at F = 0. We set up the field 
F = 0.03 au at t = 0.0. Our plots show distributions at t = 0.0 (full circles), 
t — 5.0 au (empty circles), t = 10.0 au (full squares). A total of 10^ runs 
have been performed. 

Figure 2: upper, a^x] lower cijon versus impact velocity. Full diamonds, 
F = 0.0; open circles, F = 0.03 au ; full circles, F = 0.04 au . 
Figure 3: upper, (Xcx', lower (Tion versus impact velocity. Full diamonds, 
F — 0.0; open circles, F = 0.03 au; full circles, F = —0.03 au; open squares, 
F = 0.03 au (orthogonal). 

Figure 4: upper, acx', lower aion versus impact velocity. Open circles, F — 
0.03/2"' au; full circles, F — —0.03/2'^ au. The electron is in the initial status 
n = 2. 

Figure 5: upper, charge exchange; lower, ionization. In ordinate appear 
scaled cross sections cr/n^; in abscissae scaled velocity v n. Open circles, 
F — 0.03 au, n — 1\ full circles, F — 0.03/2^ au, n — 2; open squares, 
F = —0.03 au, n = 1; full squares, F = —0.03/2^ au, n = 2. 
Figure 6: upper, d(7cx/db; lower, daion/db versus impact parameter b. Im- 
pact energy is kept fixed at E = 50 KeV/amu. Full diamonds, F = 0.0; open 
circles, F — 0.03 au; full circles, F — —0.03 au. 

Figure 7: upper, acx', lower aion versus impact velocity. Field strength is 
kept fixed at 0.03 au. Full triangles H^; full circles He^^; open circles Li^+; 
full squares C^+. 

Figure 8: upper, charge exchange; lower, ionization. In ordinate appear 
scaled cross sections a/Z]in abscissae scaled velocity v/Z'^/^. Open triangles, 
0^+; the other symbols are the same as in Figure 7. 

Figure 9: distribution of a^x over n. Upper, C^"*"; lower, O^^. Open squares 

= 50 Kev/amu; full squares = 100 keV/amu. 
Figure 10: upper, a^x] lower ai^n versus impact velocity. Full diamonds, 
F — 0.0; open circles, F — 0.03 au; full circles, F — —0.03 au; open squares, 
F — 0.03 au (orthogonal). The length L has been kept fixed at L = 10 au. 
Figure 11: distribution of acx over n for C^~^ — H impacts at an energy 
E = 100 KcV/amu. Diamonds, F = 0.0; open circles, F = 0.03 au; full 
circles, F = —0.03 au. The length L is kept fixed, L = 60 au. 



13 



